On spatial non-homogeneity in iron pnictides: formation of stripes 
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The heterogeneous coexistence of antifcrromagnetism (SDW) and superconductivity on a meso- 
scopic scale was observed in iron-pnictides in many recent experiments. We suggest and discuss the 
scenario in which the heterogeneity is caused by formation of domain walls inherent to the SDW state 
of pnictides at a proper doping or under applied pressure. Superconductivity would emerge from 
the modulated SDW structure. The phenomenon is akin to the FFLO-phase in superconductors. 



The new vast class of layered HTSC materials, iron 
pnictides manifests diverse and uncommon magnetic, 
structural and superconducting properties [2] . Below we 
address competition between magnetism (SDW) and su- 
perconductivity (SC) that take place near the so-called 
"quantum critical point" (QCP) on the pnictides' phase 
diagram tuned by either doping, or pressure. In experi- 
ments 0, 0] the competition assumes a form of the het- 
erogeneous phase coexistence. We demonstrate that the 
heterogeneity can be due to the new SDW phase, where 
the staggered magnetization is modulated by the emer- 
gent periodic stripes, the lattice of solitons. In this vicin- 
ity SC arises on the background of the Soliton phase (SP). 
Density of states (DOS) in SP has the same order of mag- 
nitude as in the paramagnetic (PM) phase. 

Currently, the consensus is that the weak coupling 
nesting model Q correctly describes the most typical pe- 
culiarities of the pnictides', at least qualitatively. 

We confine ourselves to the picture of only two Fermi- 
surfaces (FS): the one for holes at the F-point, (0,0) - 
point and the other, for electrons, at the M-point, (tt, tt), 
in the folded BZ. The model is in broad use in the liter- 
ature (see, for instance, in Q) and among other things 
reproduces the overall view of the phase diagram as func- 
tion of doping 0, the interplay and the sequence of 
magnetic and structural transitions in the undoped com- 
pounds In particular, by way of changing the degree 
of nesting the model provides the built-in mechanism for 
the competition between magnetism and superconductiv- 
ity at doping or under applied pressure. 

Among many aspects of the original nesting model [1] 
that were recently investigated afresh in numerous the- 
oretical papers, there exists one interesting feature that 
deserves more attention. The phenomenon consists in ap- 
pearance of a spatially non-uniform SDW state in pnic- 
tides under rather general conditions near QCP. 

Our attention to such a possibility was attracted by 
experiments (sl, [ll- In Q a spontaneous spatial hy- 
brid SDW/SC structure (of few nm) was reported in 
the SrFe2As2-crystal under pressure. In Q the hetero- 
geneous phase coexistence was observed in single under- 



doped (Ba,K)Fe2As2 crystal both above and below SC 
Tc- The spatial scale in Q amounted to 65 nm. 

In the scenario under investigation the tendency to 
form the periodic domain structure is inherent into the 
SDW nesting mechanism itself. Therefore, onset of the 
modulated SDW phase may occur at temperatures even 
above SC Tc- Precisely such behavior was found in the 
quasi-one-dimensional Bechgaard salt, (TMTSF)2PF6, 
The modulated SDW (at a fixed pressure, P < Per, 
i.e., before QCP) was seen above and below the transi- 
tion into the new SC state and interpreted in terms of 
the new ground state: SP 12 1. 

First, as it seems, modulated CDW or SDW phases in 
three dimensions were discussed in 13, l3|- In 14 1 the 
authors investigated the model [1] of the two anisotropic 
Fermi surfaces with the shapes deviating from the perfect 
nesting. Consequently, the energy spectrum for electron 
and holes was chosen as iiTe.ft = ivFi,2{t — ??i,2(?^)), 
correspondingly. Deviations from the ideal nesting rjii^) 
and ?72(~?^) describe both the anisotropy and the doping. 

Analysis [14] had shown that the instability line with 
respect to transition into a commensurate CDW (or 
SDW) phase in the (T, 77) -plane possesses the reen- 
trant character, similar to a sup erconductor placed into 
the exchange field, loz [1, [l^- On the side of large 
("antinesting") terms, r]i{j^) and ?72C^), the onset of 
SDW first occurs at an incommensurate (IC) vector 
where QoiT^,''^) is the commensurate vector 
at zero ?7i(~p^), r]2{l^)- At T = the g-vector is related 
to the magnitude and the anisotropy combined, as given 
by the "antinesting" factor, |77i(?^) — 'i]2C^)\ [l4j : 



(Ini I (771 -772)^ 



(cost/?) ^(7^ I 



■^Tcq(Vfi + VF2) J 

(1) 

(Average, <...>, is taken along FS). 

Appearance of an IC SDW on the phase diagram of 
pnictides was discussed earlier (See for instance [a. 
However, we emphasize the profound difference between 
IC SDW and the inhomogeneous SDW state, where elec- 
tronic degrees of freedom inside domain walls (DW) are 
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FIG. 1. The (r, 77) phase diagram {rj - tuning parameter): the 
instabihty line (dashed) for commensurate SDW shows reen- 
trance below the tri-critical point {Ttcr, Vtcr)', SP starts before 
the putative 1-st order transition line (dotted) and extends to 
the PM phase boundary (dash-dotted). The QCP is covered 
by the SC dome. See text for the details. 



not gapped which leads to a finite and spatially inhomo- 
geneous DOS (see below). 

On the (r, 77)-phase diagram in Fig.l the instability 
line (dashed) for the commensurate SDW {q = 0) dis- 
plays reentrance below the tri-critical point {Ttcr,T]tcr)', 
the dotted line that goes down from the same point is 
the would-be line of the 1st order transition between the 
commensurate SDW and the paramagnetic phase. The 
periodic SP supersedes this transition. The new phase 
extends to the dash-dotted line on the right. The posi- 
tion of the QCP inside the SC dome is shown tentatively. 

We derive the system of the Eilenberger-type equa- 
tions describing the inhomogeneous SP state. With 
the view to demonstrate the new phenomenon, we ac- 
cept the model [ij] assuming a density-density inter- 
action between the e- and /i-pockets. The analytical 
solution for SP is not available. In particular, if the 
ideal nesting is broken by doping only, the substitution, 
vprjii^) = —VFrj^il^) — const ^ I, leads to equations 
identical to the ones for FFLO problem [l^], solved 
numerically in Same methods apply to the equa- 

tions for Fe pnictides in the general case. Later on in the 
paper, we investigate emergence of domains from the ho- 
mogeneous SDW/CDW state at the increasing strength 
of the parameters, ?7i(^), ^(l^), say, by pressure. 



Derivation of the Eilenberger-style equations for the 
itinerant 3D or 2D SDW/CDW case is similar to their 
original derivation [l8| from the Gor'kov equations in su- 
perconductors and is described below only briefly. (Sim- 
ilar equations were used for ID physics in [l9l]). 



Introduce the matrix Green function 

Kj\X, X j I ^ ^ 

\ — (J21 —(-^22 



(2) 



Here Gik{x,x') — (t {if) ^{x)il}^ {x')^ , x = {r,T), and the 

non-diagonal indices, 12, 21, belong to the ofF-diagonal 
Green functions, non-zero in the presence of SDW/CDW- 
order. The two FS's are connected by the structure vec- 
tor To save on the notations, we consider the CDW 
case. The only difference is in the spin matrix, {a ■ I), 
present as a factor, in the definitions of the non-diagonal 
SDW components, Gi^kix, x') {I stands for the staggered 
magnetization direction). 

We then write down the first pair of equation: 



{d/dr + H,)Gn [x, x') + A^( f )G2i {x, x') = 1 
{d/dT + Hu)G2i {x, x') + A*Jr)Gn {x, x') = 1 



(3) 



Here: 

A^(f) = AG2i(r;ri = ra); Al^(r) = AGi2(r;ri = T2) 

(4) 

( A is proportional to the interaction constant). 

To account for the spatial dependence in ([4]), we write 
for free electrons and holes: 



e,h 



^(V /2mi,2 + ^J.) T VFi,2 ■ rjiAP 



(5) 



The "gap" paramete rs (111 ) couple the electron and the 
hole FS's. Following [18|, introduce the operator in the 
L.H.S. of Eqs.dSl): 



d/dr + H, -A^ir) 
Al^(r) -d/dr-Hh 



(6) 



so that G^^ X G(a;, x') = S{x—x'). The operator that acts 
on the matrix ^ from the right side: G{x,x') x G^^ = 
d{x — x') is obtained by changing d/dr —d/dr. 



Substitute d/dr 



On in G^^ and d/dr 



in G^^ and rewrite equations for the frequency Fourier 
components. The Green functions are not diagonal in 
the momentum representation: 

G{x,x') G{uj„;r,r') =^ Giujn;p,p- k) (7) 

The essence of the "quasiclassical" method [3, [2^ in 
the theory of superconductivity is in integrating out the 
energy variable = vpt = vf{p — Pf) thus rewriting 
equations in terms of new functions: 



^G{ujn;p,p- k) 

ITT 



g{^n;pF,k) 



5ii(w„;pF,fc) 5i2(w„;pF,fc) 
-g2i{ujn;pF,k) -g22(Wn;PF,fc) 



(8) 
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In a superconductor, the characteristic scale for Tc's 
is usuahy much less than the energy Fermi: Tc <SC Ep. 
With the characteristic scale for the magnetic transition, 
TsDW ~100-200 K and the Fermi pockets energy scale 
~ 0.2 eV [3], the quasiclassical approach is expected to 
work well enough for pnictides as well. 

The equations for the integrated Green function in 18, 
[ioj l emerge from commuting the matrix equations for the 
Green functions: G^^ x G — G x G^^ ~ from whence 
the variable drops out. Omitting the arguments of 
g{ujn ;PF,k) and denoting Sjj^ =vf [ryi ("pf ) -'72 ("pf )], the 
resulting equations acquire the form: 

-t{vF ■ V)<?ii + A^(f)52i - A*_^(f ).gi2 = 
{~i{vF ■ V) - 2iujn - Spp}gi2 + 2A^{r)gii ^ 
{+i{vF ■ V) - 2za;„ - (5p-^}<?2i + 2A*_^{r)gii = 

_ , (9) 

As in [18| . gii = (722 and the normalization condition is 
gfi — (?i252i = 1- The self-consistency Eqs.Q now is: 



A.(f)=AT 



dS^^^g2i{uJn;PF,k) (10) 



Eqs.(|9l ITOl) form the closed system. (One drawback 
of the "quasiclassical" method is that at the derivation 
one needs vfi — vf2, "mi — ni2)- In Fe pnictides two 
FS are connected by the commensurate vector ^o(7'',7r). 
Therefore A(r) is real and we omit the index Q below. 

Parameters ry depend on the specific material. There- 
fore we stay below only on the principial properties of 
the SP. Near the right boundary of the area on the phase 
diagram in Fig.l occupied by SP, the amplitude of the 
SDW- parameter. A, is small and the periodic solution 
of the form A(r) = A cosfgV) can be constructed by per- 
turbations in Eqs.([3l|4]) [e[i3]- With the decrease of the 
?7-parameters and entering the developed SP, the period 
gradually increases; below a threshold at a critical value 
of the ?7-parameters, the system enters the homogeneous 
SDW phase. Slightly above the threshold the structure 
consists of the almost isolated domains. 

It presents the decided interest to consider formation 
of one single DW at the threshold as it sheds more light 
on peculiarities of SP. Assume the dependence of A(a;) on 
one spatial variable. Return to Eqs.Q and address the 
problem of the eigen values and the two-component eigen 
functions {u,v) for the operator ([5]) [d/dr —E): 

{-ivF,xd/dx - vftii{pf))u + A{x)v = E{vf,x,Pf)u 
{+ivF,xd/dx + VFri2{PF))v + A{x)u = E{vf,x,Pf)v 

(11) 

The substitution w, u =^ (u, i;) xexp(z[ ^"^ ('?i+?72)]a;) 
transforms into: 



—ivF,xdu/ dx -\- A( x)v — E{vf,x)u 
+ivF,xdv/dx + A{x)u — E{vf,x)v 



(12) 



where E{vf,x) = E{vf,x,Pf) + Spp/2. 

Consider the energy spectrum for system (jl2p in the 
presence of a rarified periodic array of DW. Near the 
isolated wall choose the "gap", A{x), in the form: 



A(x) = Ao • tanh(a;/(a^o)) 



(13) 



The distortion to the homogeneous state (A(a::) = Aq) 
thus produced, is energetically unfavorable. At large sep- 
aration between the walls the resulting energy loss is the 
sum of contributions originated near each domain, ^ Eg . 
(Here Eg is the one wall energy cost per unit length). 

Return now to the putative 1"* order transition de- 
picted on Fig.l. To the homogeneous SDW/CDW phase 
corresponds the gain in the Free Energy density: 



AFsDW = -2iy{eF)Al 



(14) 



{i^{sf) — m/4:TT - DOS per single spin). But deforma- 
tions, r]i and 772, of the two initially congruent FS's also 
decrease the energy of the normal (PM)state: 

AF„ = -2iy{sFm^Fm)^) + {{vFm)^)] (15) 

Hence, if for instance, the system is doped, the last 
mechanism may offset the energy cost. Eg > 0, of forming 
single DW owing to the fact that there are locally gapless 
states inside the wall in the form (IT^ . Parametrizating 
doping, as before, vftJi = —vf'I]2 — /, one can write the 
local PM contribution of Eq. ([T5)) in the form: 



(16) 



AF„(x) = -AvieFpix/xo)!^ 



where xq is a scale for the domain width. Denote by 
^normi^) integrated contribution from one single do- 
main (per unit wall length): 



f: = 



-AvieF)!"^ j v{x/xo)dx (17) 



The threshold value, /c, is determined by: 







(18) 



At / > /f, the spontaneous formation of domains will be 
arrested by their interactions. 

In principle, finding Eg, vixjxQ) and the very profile 
of A(x) can be reduced to the self-consistent mean- 
field problem by solving Egs. pTj) for the band energy 
spectrum and the eigenfunctions in the periodic poten- 
tial, A{x + L) ~ A{x). Near the wall where A(a;) has 
the form (fT3l) Eqs. (fT2|) a.t L ^ xq can be solved in terms 
of the hypergeometric functions. However, the need to 
account the periodicity of (u, v) at large separations be- 
tween walls makes the approach rather cumbersome. 

With our main goal to attract attention to this new 
phenomenon we choose the model [l3| that allows the 
one-to-one mapping on the FFLO problem. Therefore it 
is possible to avoid these tiresome calculations by instead 
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making the use of the numerical results 17[. Without 
staying on further details, we obtained for v(x/(,o)) ~ 
1.2c/i-i(0.63a;/^o)- Note that V{x/^o)) ^ 1 for x of the 
order of i-e- local DOS inside domain has the same 
order of magnitude as in the PM normal phase. That 
is also true in the periodic soliton lattice. With Ic = 
0.655Ao we found from Eq.((T8]) the soliton energy 



(19) 



Es = 10.27 • HeF)CaAl = 0.26 • ppAc 



(We accepted the BCS value — fivp /t^^q- The Planck 
constant h is restored in the expressions for i^o-) 

With the soliton energy known. Eg . ([T9l) . turn now to 
the role of the anisotropy. Assume / 6{pF)dlp = 0, the 
equal numbers of electron and holes. Coming back to the 
relation between E{vf,x) and E{vf,x,Pf) in Egs. ljlll) and 
(IT2]) . one sees that the energy spectrum of Eas. p2)) for 
(m, v) is pertinent only to distorted homogeneous phase, 
when negative contributions of the form of Eq. ()16|) are 
absent. The energy spectrum of Eqs. dT^ . hence, con- 
tributes only to calculations of the soliton energy Es- 
Account of the (Jjj^-terms, E{vf,xtPx) = E(yF,x) —5ppl2, 
decreases the cost of the single wall by filling up the mo- 



mentum states with Sj^p > 0: 



Ek 



(Uf 
2^ 



(20) 



In the notations t(p||) = = ppcosLp, t(p||) now 

defines the energy spectrum for carriers moving inside 
the wall. For the corresponding DOS i^wi^) (per single 
spin, per unit length) one has: 



(21) 



i^wie) ^ / 5{e -t{p\\))dp\\/2Tr 



By the order of magnitude Vw{£) ~ {dt{p\\) / dp\\) ^ ^ 
pp/t; with t - Ao, ly^ie) ^pf/^o- The DOS in Eq.jll]) 
is concentrated inside the domain width ~ ^o- Therefore 
the local DOS in the periodic soliton lattice with L ~ 
is large, i.e., again of the same order as v{£f) = m/ATT 
the 2D DOS in PM state. 

Thus Ekin = Eg determines the threshold for onset of 
the modulated SDW state at large enough geometric mis- 
match between the two FS's. (The mechanism is akin to 
the one for forming periodic domains in (TMTSF)2PFg 

HEi.) 
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Up to this point we studied only the variation of the 
SDW state with changes to the degree of nesting. Mean- 
while, reducing the SDW gap also opens way to emer- 
gence of a SC phase. Interactions responsible for the 
1** order magnetic and structural transition seem to be 
stronger than the ones that lead to the SC pairing. In- 
deed, in the stoichiometric phictides the former occurs at 
higher temperatures. In the area of the phase diagram 
discussed so far SC would develop on the background 
of the SP. Note that when SC is still weak, its appear- 
ance may not fully remove remnant DOS that, as we 



seen above, has by order of magnitude the same value as 
in PM state. This may explain the substantial residual 
density of states towards K, revealed via nuclear spin- 
lattice relaxation rate in SC domains of SrFe2As2 Q and 
via the finite linear coefficient in the specific heat in SC 
state 



22|. 



The heterogeneous phase coexistence reported in ex- 
periments 3| seem to realize the discussed scenario. In- 
deed, in [3| the heterogeneous state first sets in below 
TsDW «70 K, i.e., weU above w32 K. While the mag- 
netic force microscope (MFM) images do not visualize 
a periodicity in some special direction, this, actually, is 
not expected. Below the 1** order transitions samples 
are twinned. In addition, unlike the strongly anisotropic 
(TMTSF)2PF6 [13, the spontaneous formation of do- 
mains in pnictides can emerge along any arbitrary di- 
rection leading to a pattern similar to the one seen in 
[1] . More recently, the scanning force microscopy (SFM) 
measurements [23j in CaFei.94Coo. 06^82 revealed the ex- 
istence of striped electronic nanostructures with dimen- 
sions ~ 4 nm. Characteristic spatial scale for the su- 
perstructure is ^0 ~ ^-^^Tt-vf /TsDW as obtained above. 
With the typical vp ^ 10*" cm/sec and Tsdw ~ 100 K 
for pnictides it leads to few nm. The /iSR measurements 
(2^ ] were capable to determine only the volume fractions 
of the coexisting SDW and SC sub-phases. The hybrid 
SDW/SC structure found in Q may correspond to the 
case when Tsdw is low and close to Tc, but we have no 
results when SDW and SC compete for existence on equal 
terms, i.e., under the SC dome near QCP in Fig.l. Such 
a competition remains the subject of the great interest 
(see, e.g., d,!!!]). 

To conclude, we have shown that emergence of a het- 
erogenic state with large local DOS is an intrinsic prop- 
erty of iron pnictides at low temperatures. Mechanisms 
of formation of a single DW near the threshold were in- 
vestigated unveiling the difference in the role of doping 
and anisotropy. Although the simple model [l3l omits 
many details concerning the interactions and the energy 
spectrum in real materials, it seems to confirm the view 
that mesoscopic phase separation observed by means of 
NMR, ^SR, MFM and SFM, in reality is nothing but the 
formation of the system of stripes. 
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